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ABSTRACT 

High accuracy spacecraft tracking requires tropospheric 
modeling which is generally scaled either by estimated or 
measured values of surface refractivity. This report sum- 
marizes the results of a worldwide surface refractivity test 
conducted in 1968 in support of the Apollo program. The 
results are directly applicable to all NASA radio tracking 
systems. 
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SUMMARY 


The troposphere is the neutral (i.e., un-ionized) portion of the atmosphere 
extending from the Earth's surface to a height of approximately 30 km. Tropo- 
spheric tracking data bias, which is elevation angle dependent, can at 5® elevation 
introduce biases up to 20 meters in range, 50 cm/sec in range rate and several 
milliradians in angle (Reference 1). At zenith the tropospheric range bias de- 
creases to 2 or 3 meters while the range rate and angle tropospheric biases 
tend to zero. The ionized portion of the upper atmosphere, the ionosphere, is 
also a significant bias source at VHF tracking frequencies (~150 MHz). However, 
ionospheric effects decrease inversely as frequency squared and at S-Band 
(~2 GHz) or higher are negligible during routine orbit computation. Ionospheric 
effects are not discussed in th’C paper. 

The frequency independent tropospheric biasing effect on range, range rate, and 
angle data is highly correlated with surface refractivity at the site of radiowave 
reception. Refractivity is simply a convenient means for specifying index of 
refraction where the relation between the two is given by: 

n = 1 + N(10“®) 

n = index of refraction 

N = refractivity 

In this report the measured time dependent variation in surface refractivity is 
presented for 12 Apollo sites during 2 months in 1967 and for 15 sites during 
8 consecutive months in 1968. This information is useful in estimating the ab- 
solute tracking accuracy attainable if, as generally is the case, average values 
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of surface refractivity are used in the reduction of tracking data from any given 
site, The original purpose of this test was to assure that the estimated monthly 
surface refractivity values used by the Manned Spacecraft Center (NASA-Houston) 
were adequate for the Apollo program. It was seen that Apollo tropospheric 
corrections were consistent with the mission required spacecraft tracking ac- 
curacy. All Apollo stations were assigned monthly average values of refractivity 
which in no case differed by more than 10% from the measured monthly averages. 

However, for tracking experiments such as Very Long Baseline Interferometry 
(References 2 and 3) the tropospheric variations throughout the day or month 
cannot be neglected. This is also true for the stringent tracking accuracy as- 
sociated with geodesy and Earth physics studies (Reference 4) such as now being 
undertaken at the Goddard Space Flight Center. 

The accuracy of the indivi -*ual refractivity values calculated in this report is 
estimated to be on the order of 1% (Reference 5). The resolution of these 
measurements is ± 0.3 N units or typically ± 0.1%. 
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SURFACE REFRACTIVITY MEASUREMENTS 


AT NASA SPACECRAFT TRACKING SITES 
1.0 INTRODUCTION 

It is well known that the troposphere biases spacecraft tracking data and conse- 
quently some type of tropospheric correction is required prior to orbit compu- 
tation. The best correction would of course be obtained from a measured profile 
along the ray path between tracking site and spacecraft. However, such informa- 
tion is not usually available and instead some type of mathematical model is 
assumed. Fortunately the troposphere is "well behaved" in the sense that a 
reasonable integrated profile can be computed with the only input parameter 
being surface refractivity. Reference 6, which presents a comparison of various 
tropospheric data correction equations, clearly shows the scaling effect of 
surface refractivity. For purposes of radio refraction and time delay correction 
an exponential troposphere such as presented by the National Bureau of Standards 
offers good agreement with measured profiles for most tracking applications 
(References 1 and 7). 

It should be noted that for short baseline interferometers such as the NASA 
Minitrack system there is an inherent first order correction for tropospheric 
biasing effects. This inherent correction results whenever one uses the value of 
freespace wavelength of the arriving radio signal to calculate the interferometer 
angle (Reference 8). 

The measured data presented in this report permit an estimate of the tropospheric 
errors due to variations in surface refractivity alone. There are of course other 
modeling errors which may predominate, especially at elevation angles below 
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about 10“. Horizontal gradients are not included in the usual model where a 
spherically stratified (i.e., horizontally homogeneous) troposphere is assumed. 
Also at low elevation angles tropospheric turbulence can introduce noise into 
the tracking measurements. However, these noise like errors tend to be aver- 
aged out by the data smoothing and orbit computation process and are therefore 
less serious than steady state or biasing errors. It should be noted that the 
noise considered in the foregoing sense is not a radio noise such as attributed 
to thermal emission by the atmosphere but rather a phase modulation of radio 
signals traversing regions of time varying permittivity and hence time varying 
radio index of refraction. 
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2.0 CALCULATION OF SURFACE REFRACTIVITY 


The correlation between surface refractivlty and total tropospheric biasing 
effects on tracking data is well established (Reference 6). This report is con- 
ceriicd only with presenting the means for calculating surface refractivlty based 
on surface meteorological data and showing how this value of refractivlty varies 
from site to site throughout the year. The relationship between meteorological 
data and refractivlty is also well known. The relationship reported by Smith and 
Weintraub in 1953 (Reference 9) is still in widespread use by radio scientists 
throughout the world (Reference 5). Reference 5 presents an excellent review of 
the physics relating the radio dielectric constant (and therefore index of re- 
fraction) to the parameters temperature, total air pressure, and relative humidity 
and the derivations will not be repeated here. 

2.1 SCOPE OF TEST 

Meteorological data was sent in from the 15 Manned Space Flight Network 
stations indicated in Figure 1. Station personnel measured the air temperature, 
wet bulb temperature (using various types of psychrometer) and air pressure, 
at an average time interval of every 6 hours. The test (Reference 10) was con- 
ducted throughout August and September 1967 and then again on a continuous 
basis from February 1968 to September 1968. The raw data was mailed to the 
Goddard Space Flight Center from each station and entered onto computer cards 
which were then used as input to a computer program which performed tht' p- 
propriate conversion of units, refractivlty calculations, and plotting. The results 
were excellent and based on the reasonable assumption that readings were made 
to within ± 0.25®C in temperature, the absolute error in N^ is less than 3 N units 
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(Reference 14) 

igure 1. MSFN Sites Which Participated in the Surface Refractivity Test 



(Reference 6). The corresponding relative changes (i.e., resolution) is on the 
order of 0.3 N units. The accuracy of the individual refractivity values cal- 
culated in this report is thus on the order of 1% (since is on the order of 
300). The corresponding resolution is 0.1%. 

2.2 MATHEMATICAL FORMULATION FOR REFRACTIVITY 

All data mailed in from the tracking sites included the time associated with each 
measurement. Since several different types of psychrometers, thermometers, 
and barometers were used by the various stations some dat ^ conversion was re- 
quired. Where required the following relationships were used. 

^(mb) ' ^/INCHES \ 33.86395 

i MERCURY/ 

'^(°K) = T^oc) + 273.16 

The tropospheric surface index of refraction is greater than unity and can be 
written as: 


n = 1 + N^(10"®) 


( 1 ) 


where: 

n = Surface index of refraction 
Ng = Surface refractivity 

The value of will range from a minimum of 200 to a maximum of 450 de- 
pending on geographical location, time of year, weather conditions; and so on 
(Reference 11). 
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Surface refractivlty, N,, is a function of atmospheric temperature, pressure, 
and water vapor r-^ntent (Reference 5) . 


That is; 


N 


s 



48W£ 

T 


( 2 ) 


where: 

T ^ Air Temperature or "dry bulb" temperature (T<elvin) 

P - Total air pressure (mbar) 
e = Partial water vapor pressure (mbar) 

The partial water vapor pressure is in turn related to psychrometric measure- 
ments by; (Reference 5) 

" e* - 0.00067 (T - T')P (3) 


where; 

T and P are as defined above and 

T' = The "wet bulb" temperature (’Kelvin) 

e^ = The saturation vapor pressure at T' (mbar) 

The value of e^ can be obtained directly from appropriate psychrometric charts 
(Reference 12). However, in order to facilitate computer computation, an em- 
pirical foxmulation was used which resulted in negligible error relative to the 
uncertainty in the meteorological field measurements. The expression used to 
calculate the saturated vapor pressure in millibars is given by (Reference 13) 

a t 

e; = 6.11 X iQb + t (4) 
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where; 


a = 7.5 
b = 237.3 

t = wet bulb temperature (®C) 
note that 

T* ^ t + 273. 16 (^Kelvin) 


Thus the value of N^, which can be directly related to radiowave tropospheric 
refraction and time delay, is determined by three basic measurements — namely; 

1. Dry bulb temperature; 

2. Wet bulb temperature; 

3. Total atmospheric pressure. 


Once all the data from a given station for a particular month was processed, the 
monthly mean and standard deviation were computed usfng; 


mean “ “ / — m = number of observations 

s La m 

i= 1 


standard deviation 


t. 

i = 1 

m - 1 


( 5 ) 


( 6 ) 


As a matter of interest the computer processed some 50,000 measurements and 
provided point to point 35mm film plotting as well as printout of over 15,000 
values of refractivity. The next section of this report presents the surface 
refr activity measurement results. 
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3*0 SU RFACE R EFRACT IVITY MEASUREMENT RESULTS 
3.1 MONTHLY VARIATIONS 

Table 1 summarizes the maximum, minimum and monthly average of refractivity 
as obtained from the 15 Apollo sites. The Antigua (ANG) station did not furnish 
pressure data so a nominal 76.2 cm of mercury atmospheric pressure was used 
to obtain the ANG values in Table 1. This value of pressure corresponded to 
the average at Carnarvon, Australia which is approximately the same height 
above sea level as Antigua (see Figure 1). The maximum absolute variation in 
pressure at Carnarvon throughout the test was 1.5 cm of mercury. It is there- 
fore estimated that the ANG values are correct to at least ± 5 N units. Also it 
was noted (after the fact) that the GDS August- September 1967 pressure data was 
in error. Therefore these 2 months in Table 1 were computed using P = 910 mb 
based on 1968 GDS data. 

Figure 2 shows the measured total atmospheric average pressure during mid- 
May 1968. The expected correlation between station altitude and pressure 
(Figure 2) is quite apparent. In this test the pressure values were seen to be 
primarily a function of station altitude with maximum variation of only a few 
percent of total pressure throughout the test. 

The computer plots of refractivity measurements versus time from all sites 
(except ANG) then make up the bulk of this report. These results clearly show 
the degree of variability in surface refractivity one can expect throughout any 
given month. Depending on geographical location, the one-sigma spread is seen 
to be typically ± 10 to i 15 N units over a given month. 
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Table 1 

Monthly Values of Refr activity 
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Figure 2. Example of Measured Atmospheric Pressure Distribution 


This standard deviation (labeled SD in Table 1 and the refractivity plots) is a 
measure of the variability throughout a given month for a particular station. It 
is not a measure of measurement accuracy. The maximum and minimum values 
of refractivity reflect temperature extremes, low temperatures resulting in high 
refractivity and vice versa. The plots indicate the diurnal variations which again 
are primarily temperature dependent, lower temperatures generally being as- 
sociated with local night. 
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3.2 MONTHLY AVERAGES COMPARED WITH APOLLO VALUES 


Table 2 lists the monthly average refr activity estimates used by the Apollo 
program. Table 3 then presents the differences between measured monthly 
averages and the corresponding Apollo estimates. In general the agreement is 
good (within 13% or better for all cases) considering the fact that measured 
monthly averages ranged from 266 to 385, with individual values ranging from 
231 to 434. Where Apollo refr activity averages are consistently above or below 
observed monthly averages (for example BDA or CRO) some adjustment of the 
Apollo values may be warranted. 


Table 2 

Monthly Average Refractivity Estimates Used in Apollo Program 


Apollo Monthly Averages of Refractivity 


Station 

Jan 

Feb 

Mar 

April 

May 

June 

July 

Aug 

Sept 

Oct 

Nov 

Dec 

ACN 

361 

364 

366 

363 

361 

356 

353 

353 

354 

356 

355 

357 

ANG 

360 

359 

359 

361 

372 

375 

378 

379 

376 

377 

372 

366 

BDA 

335 

333 

335 

343 

353 

372 

377 

378 

371 

359 

345 

338 

CRO 

338 

339 

341 

348 

362 

378 

385 

386 

381 

361 

349 

339 

CYI 

333 

336 

336 

335 

338 

346 

343 

352 

354 

350 

343 

336 

GBM 

343 

354 

352 

361 

370 

375 

386 

383 

384 

373 

364 

356 

GDS 

279 

281 

278 

279 

297 

277 

279 

279 

278 

277 

276 

279 

GWM 

372 

377 

378 
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381 

377 

373 

379 

381 

380 

375 

379 

GYM 

318 

322 

320 

321 

331 

351 

368 

370 

372 
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328 

321 

HAW 

299 

298 

297 

299 

301 

305 

308 

310 

309 

309 

304 

302 

HSK 

306 

312 

309 

293 

300 

297 

296 

295 

298 

301 

300 

301 

MAD 

295 

293 

295 

294 

298 

302 

299 

300 

306 

304 

299 

296 

MIL 

338 

339 

341 

348 

362 

378 

385 

386 

381 

361 

349 

339 

TAN 

310 

310 

313 

305 

298 

292 

291 

290 

291 

293 

303 

310 

TEX 

328 

333 

340 

350 

362 

377 

395 

375 

372 

353 

336 

334 
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Table 3 

Observed Minus Apollo Values of Ng 


Station 




1968 




1967 

Feb 

]Mar 

Apr 

May 

Jun 

Jul 

Aug 

Sept 

Aug 

Sept 

ACN 

-15 

-15 

-14 

-14 

-13 

-16 

-19 

-17 


-10 

ANG 


- 8 

+ 1 

0 


- 3 

d 1 

H- 2 



BDA 


- 6 


- 7 

- 7 

- 4 

- 5 

- 9 


-10 

CRO 

+19 

+ 13 

+19 

+11 

+14 

+13 

+15 

+21 

+ 9 

d 9 

CYI 

- 8 

-13 

- 8 

- 7 

- 9 

- 5 

- 9 

-12 

-18 

-15 

GBM 


-14 

- G 

- 1 

+ 9 

- 2 

+ 2 

- 2 

- 4 

- 3 

GDS 

+13 


+ 5 

-10 

+ 8 

+ 3 

+ 1 

- 8 

+ 4 

+ 8 

GWM 

- 9 


- 5 

+ 2 

+ 4 

+12 

+ 6 

+ 3 

+ 7 

+ 4 

GYM 

1 

i 7 

+ 2 

+ 6 

+ 9 

+17 

+14 

+14 

+ 6 

+11 

+ 1 

f 

1 HAW 

» 

+ 8 

+ 9 


+ 9 

d 6 

+ 8 

+ 3 

+ 6 

+ 5 

+ 5 

! IISK 

1 


... . .1 

- 2 

-12 

-17 

-17 

-17 

-20 


t 

MAD 

ft 

- 3 

+ 1 

- 9 

-15 

-14 

-34 

-13 

-17 

-20 

MIL 

-24 

-15 

- 6 

-10 

- 8 

-12 

-14 

-16 

■ 

- 9 

TAN 

- 2 

- 3 

- 3 

- 3 

- 3 

-1- 3 

+ 1 

- 1 

H 

0 

TEX 

+ 5 

+15 

+25 

+19 

+ 8 


+ 8 

+ 4 


0 


Note: 

1. Apollo values are NASA-Manned Spacecraft Center best estimates of monthly averages at 
each site. The same values are used each year. 

2. Range of measured monthly average AN^ = 119 

where maximum - 385 (TEX June 1968) 

minimum = 266 (MAD August 1968) 
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Finally, Table 4 compares the August 196V - August 1968 and September 1967 - 
September 1968 monthly average refractivity. The year to year variation in 
monthly average refractivity is in all cases seen to be less than 6% and in 
several instances less than O.SVr. 


Table 4 

Year to Year Variation in Monthly Average of Ng 
(August, September 1907-1968) 


Staiion 

■ ' j 

August 

1 9 ■ ) i 

lllgUSt 


September 

1907 

September 

1908 

\ i 
1 

ACN 


3.‘51 


o u 

• J t ^ O 

.‘'.37 

1 1 

AKG 

\'o Data 

380 

— 

No Data 

378 

i 

1 

liDA 

.‘307 

373 

6 

301 

302 

1 

1 

CHO 

3J31 

338 

7 

334 

346 

12 

CYI 

334 

343 

9 

339 

342 

3 

GBM 

379 

385 

6 

381 

382 

1 

GDS 

283 

280 

3 

286 

270 

16 

GWM 

386 

385 

1 

385 

384 

1 

GYM 

381 

384 

3 

373 

378 

5 

HAW 

315 

313 

2 

314 

315 

1 

HSK 

No Data 

278 

— 

No Data 

278 


MAD 

283 

266 

17 

286 

293 

7 

MIL 

379 

372 

7 

372 

365 

7 

TAN 

287 

291 

4 

291 

290 

1 

TEX 

376 

383 

7 

372 

376 

4 
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4.0 CONCLUSIONS 


The results of this worldwide refr activity test have shown that: 

1. The monthly average refractivity values used in the Apollo program 
are in all cases within 13% of the measured values. This is satisfactory 
in terms of required Apollo tracking data accuracy. 

2. The day-to-day variations in refractivity relative to the monthly average 
were observed to be as large as ± 50 N units. At low elevation angles 
(5® or less) this results in residual tropospheric bias errors on the 
order of several meters in range and 5 cm/sec or less in range rate 
depending on spacecriift geometry. At 45 ' elevation these bias errors 
would be reduced to 0.5 meters in range and 0.05 cm/sec maximum in 
range rate. 

3. For routine orbit computation the use of estimated monthly average 
values of refractivity for each site is satisfactory. For high accuracy 
orbit and trajectory computation in conjunction with geodesy and Earth 
physics, tropospheric measurements should be made during the time of 
spacecraft trace rg. 

4. At any given site humidity and temperature variations are the primary 
cause of changes in refractivity. Air pressure changes at a given site 
(in the absence of temperature- humidity variations) would result in a 
maximum change of only ± 5 N units per year. 

5. This test has shown the variability associated with surface refractivity, 
a scaling quantity in most mathematical models for the troposphere. 
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A necessary continuing study is the comparison of mathematical models of the 
troposphere with measured profiles (e.g,, Reference 15). 
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